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SIGNIFICANCE  AND  EXPLANATION 

Signal  processing  has  wide  application  in  engineering  and  science.  An 
important  mathematical  aspect  of  this  field  is  the  reconstruction  of  a function  from 
its  values  at  the  integers.  This  interpolation  problem  may  be  solved  in  various 
ways  e.g.,  Whittaker  cardinal  series,  or  moving  averages,  as  with  local  polynomial 
interpolation,  or  with  spline  functions  with  an  infinite  number  of  knots. 

In  this  report  we  enlarge  the  class  of  spline  function  methods  which  may 
be  used  in  interpolation.  These  splines  are  piecewise  smooth  solutions  of  some 
prescribed  constant  coefficient  differential  operator.  In  this  report  we  allow  for 
periodic  solutions  to  the  differential  operator,  that  is,  for  complex  zeros  in  the 
associated  characteristic  polynomial. 

In  addition  to  their  usefulness  in  interpolating  a function  at  all  the  in- 
tegers, our  results  also  apply  to  spline  interpolation  on  the  circle.  Our  appli- 
cation in  this  regard  extends  and  unifies  previously  known  results. 


The  responsbility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC,  and  not  with  the  authors  of  this  report. 


SPLINE  FUNCTIONS  ON  THE  CIRCLE:  CARDINAL 
L-SPLINES  REVISITED 


Charles  A.  Micchelli  and  A.  Sharma 
(Dedicated  to  Professor  I.  J.  Schoenberg  on  his  75th  birthday.) 

1.  Introduction.  Although  the  literature  on  splines  has  grown  vastly  during  the  last  decade 
[11] , the  study  of  polynomial  splines  on  the  circle  seems  to  have  suffered  neglect.  The  first 
to  study  the  subject  in  depth  seem  to  be  Ahlberg,  Nilson  and  Walsh  [1] . Almost  at  the  same 
time  I.  J.  Schoenberg  [8]  studied  the  problem  of  interpolation  at  the  roots  of  unity  by 
splines  and  its  relation  to  quadrature  on  the  circle.  For  discrete  polynomial  splines  on  the 
circle  we  refer  to  [5].  M.  Golomb  [3]  also  considers  interpolation  by  a class  of  "spline" 
functions  in  the  complex  plane  but  his  point  of  view  is  based  on  minimum  norm  properties  of 
spline  functions.  Perhaps  the  reason  for  this  neglect  may  be  attributed  to  the  fact  that  one 
can  pass  from  the  circle  to  the  line  by  means  of  the  transformation  z -*■  exp  2tix.  This  changes 
the  problem  on  the  circle  into  periodic  interpolation  on  the  line  with  the  difference  that 
instead  of  interpolation  by  piecewise  polynomial,  we  now  consider  piecewise  exponential  poly- 
nomials with  complex  exponents. 

Recently  J.  Tzimbalario  [13]  has  brought  out  the  close  affinity  of  the  problem  of  cardinal 
trigonometric  interpolation  and  that  of  interpolating  at  the  roots  of  unity  on  the  circle  and 
has  given  a unified  treatment  using  the  ideas  of  cardinal  L-splines  due  to  C.  Micchelli  [6] . 

The  object  of  this  note  is  to  bring  this  affinity  into  better  perspective  by  considering 

n ^ 

piecewise  functions  of  the  form  £ c z which  is  suggestive  of  the  Muntz  approximation 

0 V 

theorem  (D.  J.  Newman  [7]).  Such  polynomials  have  been  termed  "incomplete"  polynomials  by 
G.  G.  Lorentz.  In  this  case  we  are  led  to  investigate  cardinal  L-splines  with  complex  exponents 
occuring  in  conjugate  pairs.  Our  result  in  this  context  extends  some  of  the  results  of  [6] 
which  treats  L-splines  with  real  exponents. 

The  possibility  of  such  an  extension  was  first  envisaged  by  Schoenberg  ([9]  p.  274)  and 
later  investigated  by  Tzimbalario  and  sharma  [12]  for  the  case  of  trigono,..jtric  splines.  Our 


treatment  is  based  upon  ideas  from  [6]  and  19] . It  does  not  require,  as  in  [13] , a lemma  of 
Ahlberg,  Nilson  and  Walsh  concerning  the  sign  of  Re  p(z)  on  the  locus  Im  p(z)  = 0 for  a 
polynomial  with  real  negative  zeros.  (See  [13]  for  its  precise  formulation). 

In  Section  2 we  formulate  the  interpolation  problem  on  the  circle  and  state  our  main  re- 
sult. In  Section  3 we  provide  a necessary  and  sufficient  condition  for  the  solution  of 
the  interpolation  problem.  Section  4 is  devoted  to  some  formulae  from  cardinal  L-splines  and 
their  relationship  to  interpolation  on  the  circle.  Here  we  formulate  in  two  different 
ways  the  uniqueness  criterion  of  Section  3.  Section  5 deals  briefly  with  the  deriva- 
tion of  the  B-spline  representation.  In  Section  6 we  obtain  results  about  the  zeros  of 
cardinal  L-spline  with  complex  exponents  thereby  extending  some  of  the  results  of  [6] . Section 
7 deals  with  a proof  of  Theorem  1 of  Section  2.  We  turn  to  an  application  of  these  results  to 
quadrature  in  Section  8 where  we  also  determine  the  generalized  monospline  of  least  L^-norm. 

2.  Statement  of  the  Problem  and  the  main  result: 

Let  U denote  the  unit  circle  |z|  = 1 and  let  n,  k be  positive  integers  k > n >_  1. 

Let  A = {X.  < X,  <...<  X } be  a set  of  n + 1 integers  which  is  a subset  of  {0,  l,...,k-l). 
0 1 n 

We  shall  denote  by  it.  the  class  of  polynomials  P(z)  given  by 


P(z)  = l 


Let  u)  = 1 (a  = ■—)  be  the  k*”*1  root  of  unity.  We  shall  consider  the  class  5-^  ^(=  S)  of 

A-splines  S(z)  defined  by  the  following  two  conditions: 


(i)  S ( z)  r Cn_1(U) 

(ii)  S(z)  | = Py(z)  e tt^,  v = 0,  l,...,k-l 

v 

where  A denotes  the  arc  of  the  circle  (u>v,  uiV+^)  . We  can  now  formulate  the 

v 

Problem.  Let  ij i = e1<T,  0 ^ t < 1 be  a point  on  the  circle.  For  what  choice  of  k,  A,  T 
does  the  interpolation  problem  ( I .p) : 

(2.2)  St^w3)  = w.,  j = 0,  1 , . . . ,k-l 


have  a unique  solution  3 


* ? 
A ,n 


■i  I 


i v fvV-  ■ ■ i 


I.  T = 0 and  n is  an  even  integer 


II.  T = — and  n is  an  odd  integer  . 

We  get  more  information  from 

Theorem  2.  In  each  of  the  cases  of  non-uniqueness  in  Theorem  1,  the  I.P  admits  a solution 
S(z)  e S if  and  only  if  w^  satisfy  the  relation 

k-1 


(2.3) 


y u-hj  w.  = 0,  h = (k  + A + A )/2  . 

.*•„  1 On 

3=0  J 


3.  A Uniqueness  criterion  for  I.P. 

We  begin  with 

Lemma  1.  For  every  integer  r,  0 <_  r k-1  there  exists  a non-trivial  solution  S^fz)  e S 
(unique  up  to  a multiplicative  constant)  which  satisfies 


(3.1) 


s^tuz)  = w Sr(z)  , z e U 


The  functions  SQ,...,Sk_1  form  a basis  for  S. 


Proof.  It  is  easily  seen  that  the  conditions 
(3.2) 


a)1*  P ^’(u)  = ur  p^’d),  j = 0,  1, . . . ,n-l 


determine  up  to  a constant  a unique  P « ir  When  r ^ A,  we  may  choose  the  normalization 


(3.3) 


n „(n)  . , r (n) 

(1)  P (ou)  - ui  P (1)  = 1 


and  for  r « A,  P(z)  = zr  satisfies  (3.2).  Clearly  p may  be  extended  from  the  arc  (1,  u) 
to  the  unit  circle  by  means  of  (3.1)  and  the  resulting  extension  S^(z)  is  in  S.  The  linear 
independence  of  Sg,...S  ^ follows  directly  from  (3.1).  We  now  show  that  these  functions 

span  5. 
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S.  Then 


Since  S (wz)  = u>  S(z),  the  first  half  of  Lemma  1 implies  the  existence  of  constants  d 


such  that  S (z)  = d S (z) . This  proves  that  S , S 


form  a basis  for  S 


Lemma  2.  The  I.P  has  a unique  solution  in  S if  and  only  if  S 


Then 


for  v e A,  the  lemma  is  proved 


then 


Cardinal  L-splines  are  related  to  a given  differential  equation 


where  T denotes  the  zero  set  of 


L-splines  are  piecewise  exponential  polynomials,  i.e 


a solution  to  (4.1)  but  globally  in 


e be  the  unique  element  of  the  form 


which  satisfies 


We  shall  write  A (x; X)  for  brevity  when  there  is  no  ambiguity.  A (x ; A)  is  a rational  func 


tion  in  A and  is  an  exponential  polynomial  in  x.  We  list  below  some  of  the  known  properties 


of  A (x»A).  (See  [6]  where  they  are  given  for  T real,  but  these  algebraic  properties  hold 


also  for  T complex) 


(a)  The  extension  of  A (x;A)  from  (0,11  to  » by  the  ecuation 


mSSMSm 


(b)  On  [0,1], 

(4.4) 


xt 


VX;A>  = [t0,  tir...,tn)  — 


e -X 


where  [t  , ...t  ]f(t)  denotes  the  divided  difference  of  f at  the  nodes  t , ...t  with  re- 
0 n u n 

spect  to  the  parameter  t. 

(c)  If  A = | A | e1U,  -u  < u <_  it,  6^  = In | A | + i(u  + 2n £)  ,1  » 0,  ± 1,  ± 2, . . . , then 


(4.5) 


o»  (x-l)6.  , ® 2ti!lx 

. r e l ,x-l  r e 

An‘X;X)  =Ip„+1<V  = 1 


- wv 


(d)  If  the  zeros  of  p , (x)  are  symmetric  about  the  origin,  i.e.,  T = -T,  then 
n+1 


(4.6) 

If  we  set 

(4.7) 


A (1-x; A |t)  = (-l)n_1A_1A  (x;A|t)  . 

n n 


n t 

n ( A ;x  1 T)  = n (A  ;x)  = r ( A)  A (x;A|t),  r(A)  = Il(e  V-A) 
n n n 0 


then  n (A;x)  is  a polynomial  of  degree  n in  A,  if  x e (0,1)  and  of  degree  n-1  if 


x = 0. 


I T 

We  now  choose  T = iaA  where  A = (A„,...,A  }.  Then  for  A 4 e , 

0 n 


xt 


(4.B) 


An( x;X)  = liaXQ, . . . ,iaXnl  — 


e -X 


Choose  X = e = w , r J A so  that  in  (4.5)  6p  - i(ra  + 2tt£)  , and  set 
(4.9) 


- # r ) , . „ , r . iax 

A (z?u>  A)  = A (x; a)  ),  z = e 
n n 


Then  from  (4.5)  we  have  for  z £ A 


(4.10) 


“ ks+r 

A (zswr|A)  = u r(ia)  n \ 

n _a>  n 

n(ks+r-A  ) 
v=-0  v 


and  because  of  (4.3),  A (z;ui  | A ) > S and  satisfies 

n 


(4.11) 


u 


According  to  Lenina  1 we  have 


A (lDZ;our|A)  = ur  A (z:ti>r  I A) 
n n 


S (z)  = A (z,-ur|A) 
r n 


When  A = {0,  l,...,n-l},  has  been  called  an  r-f lower  by  Schoenberg  [8]  in  view  of  its 
rotational  symmetry.  In  this  case  on  using  (4.8)  with  X = u>  and  taking  r = n+1  we  get 
except  for  a non-zero  constant  factor  a simpler  derivation  of  formula  (6.17)  in  [8]  which 
gives  the  polynomial  component  of  an  r-flower  on  the  arc  A^. 

Using  the  above  formulae  we  now  reformulate  Lemma  2 in  a more  useful  form. 

Lenma  3.  The  I.P  is  solvable  if  and  only  if  one  of  the  following  (equivalent)  conditions  is 


satisfied: 


n (ur;x|iaA)  4=  0,  r 4 A , t e [0,1) 
n 


N = l 


+ 0,  r j A,  1 1/ 


n (ks+r-X^) 
v=0 


5.  B-splines  for  S. 

Condition  (4.12)  above  is  reminiscent  of  Schoenberg's  approach  in  [8]  for 
A = {0,  l,...,n-l)  and  may  be  derived  in  a similar  fashion  from  the  B-spline  representation 
for  S(z)  e 5.  To  this  end  we  define  the  polynomial  VL  (z)  € TI,  by 


'fi.  (z)  = 1 

A 


X,  X 

1 n 

z . . . . z 

1 1 

X,  X 

1 n 


* V[V\ Vl1 


where  V[X„, ,X  ,1  denotes  the  Vandermondi an . (z)  has  the  property  that  the  coefficient 

. 0 n-1  A 

A , . 

of  z n is  unity  and  that  (1)  =0,  v = 0,  1, ,n-l.  With  these  polynomials  we  can 

easily  form  the  B-splines  M. (z)  for  the  class  S.  We  can  also  show  as  in  [8)  that  M,(z) 

A A 

n 

has  the  support  U A.  and  that  every  S(z)  e S can  be  written  uniouely  in  the  form 
i=0  x 


Stz)  = £o"lavVz“ 


where  on  the  arc  A . we  have 


where  the  coefficients  d are  given  by  the  identity 


When  n - 1 , it  is  easy  to  see  that  A 


Also  in  this  case 


If  we  define  convolution  f * g by 


then  we  can  see  that  M, (r)  = M 


The  Fourier  series  for  M.(z)  is  given  by 


L.  S'; 


k-1 


M^fi/iw-’)  = I aiv:*,  j = 0,  l,...,k-l 


v=0 


so  that 


A 1 k_1  ■ _ • 

S = k ^ ' v = 0,  l,...,k-l  . 


j=0 


Using  (5.1)  we  have 


” K I K 


'v  “ v+ks’ 

S=— “ 


a non  zero  constant. 


For  v € A say  v = A , we  have  from  (5.2) 


el 


C,  = * r b,  +0  . 


For  v 4 A,  we  get  from  (5.2) , 


A n X£~v 

C - K n (1  - u ) N 
v . _ v 

£=0 


where 


(5.3) 


OD 

l * 

L n 


ks 


i!  (v+ks-A  ) 
£-0 


Following  the  reasoning  in  Schoenberg  [8J  we  would  be  led  to  another  derivation  of  (4.12)  of 
Lenina  3. 

A direct  examination  of  (4.13)  of  Lemma  3 in  the  general  case  seems  to  be  difficult.  In 
the  following  Section  we  prepare  the  ground  for  proving  Theorem  1 by  using  (4.12) . 


6.  Zeros  of  II  (A;x|t)  , T = {t„,  t,  , ,t  }. 

n 0 1 n 


If  T is  real  the  properties  of  the  zeros  of  n (X;x)  are  known  [6],  When  T is  com- 

n 

n 


plex  but  pn+1(x)  = IHx-t  ) is  real,  we  can  still  prove 
" 0 v 


Lenina  4.  If  the  polynomial  Pn+1<x>  is  real  with  |lm  t^|  < ir,  then  for  any  A < 0,  the 
exponential  polynomial  ( X ; x)  has  exactly  one  simple  zero  in  [0,1). 

Proof . Observe  that  if  (D-a) (D-a)g  = 0 with  |lm  a|  < n,  then  g(x)  can  have  at  most  one 
zero  in  [0,1] . 
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. Ami  . . 


Assuming  that  pn+1 (x)  has  a complex  zero 


a 


< ii 


it  follows  from  the  fact  that 


(D)  - (D-o)  (D-a)  h(D) 


with 


The  existence  of  a decomposition  as  in  (6.2)  follows  from  the  identities  (D-t)y  ■ e 


Now  suppose  to  the  contrary  that  n (X>x)  has  more  than  one  zero  in  (0,1)  counting  multi' 


plicities.  Then  according  to  (4.3)  and  (4.7),  A (x,*X)  has  at  least  2n  zeros  in  [0,n)  . Thus 


g(x)  = h(D)  A (x; X)  has  at  least  2n  - (n-1)  « n + 1 zeros  in  [0,n),  whence  it  follows  that 


there  is  at  least  one  interval 


h(D)  is  a differential  operator  with  constant  coefficient  and  A (x; X)  satisfies  (4.3),  the 


function  g(x)  has  at  least  two  zeros  in  10,1).  However  (6.1)  implies  that  (D-a)(D-a)g  - 0 


and  our  initial  remark  says  that  g(x)  can  have  at  most  one  zero  in  [0,1) . This  contradic- 


(x)  has  only  real  zeros,  we  may  proceed  as  before  and  use  the  fact  that  if 


(D-a.) (D-a,)g  =0  (a.,  a,  real),  then  g(x)  has  at  most  one  zero  in  ( 


proved  by  Micchelli  ([6)  Corollary  2.2,  p.  210) 


Theorem  3.  If  the  polynomial  p . (x)  » n(x-t. ) 


riven 


and  negative 


Moreover  these  zeros  are  simple  and  increasing  functions  of  x e (0,1) , i.e 


< 0 are  the  zeros  of  D (X;0)  then  X. (0) 


-9- 


Remark . If  all  the  zeros  are  real.  Theorem  3 is  proved  in  [6]  p.  222.  It  is  also  shown 

there  that  when  T is  real,  the  zeros  of  An(x;X|T)  interlace  the  zeros  of 

A . (x ; A I T - {t  }) . This  fact  leads  to  the  useful  information  that  each  A . (x)  is  a strictly 
n- 1 n 3 

decreasing  function  of  each  t^.  For  the  special  case  T = {0,  t,  2t,...,nt},  t real  the 
monotonicity  of  A_.  as  a function  of  t has  interesting  ramifications.  It  is  related  to 
the  problem  of  finding  best  bounds  on  the  local  mesh  ratio  of  knots  which  will  allow  interpola- 
tion by  bounded  splines  to  bounded  data  on  the  real  line  ([2],  [6]).  No  such  results  are 
known  when  T is  complex  with  T = T. 

When  T is  real,  another  derivation  of  Theorem  3 is  given  by  Schoenberg  [9]  which  is 
based  upon  the  total  positivity  properties  of  B-splines  for  cardinal  L-splines  and  upon  results 
regarding  the  generating  function  of  Polya  frequency  functions.  This  method  can  be  adapted  to 
the  case  when  T is  complex  and  T = T as  is  shown  below . 

The  proof  of  Theorem  3 follows  the  approach  of  Schoenberg  [9]  almost  verbatim.  The 
difference  rests  on  Lemma  4. 

Proof . We  shall  use  the  following  representation  for  ( X ; x)  (16]  p.  222). 

n t 

(6.4)  IT  <A;x)  = (-1) n l Q (v+l-x)AV 


where  Q , (x)  is  the  forward  B-spline  with  respect  to  the  operator 
n+l 


(-1)  Pn+1(~D)  = IUD  + t.)  . 

Since  p ,,(-x)  is  also  real  and  I Im  (-t.)l  < it,  p ( — D)  has  a Polya  factorization.  More 
n+l  j n+l 


precisely 


p Al(-D)f  = w D(w  *w  .D)  •••  (W 11W  D)  (w^f)  , w.  >0, 
n+l  n n n-i  l u u 3 J 


From  a theorem  of  Karlin  ([4]  p.  527,  Theorem  4.1),  the  sequence  {Qn+1<v+l-x) }Q  is  a Pclya 
frequency  sequence  for  all  x £ (0,1).  It  follows  from  a theorem  of  Schoenberg  (Theorem  5.3 
p.  412  [4])  that  rtn ( X ; x)  has  all  its  roots  real  and  negative. 


: ) 


, k > • A-r  ‘ V iH'fili'HSiiiir  ill' 


'Ji 


l 


Inequality  (6.3)  (the  simplicity  of  these  zeros)  follows  from  Lemma  4 on  repeating  the 
reasoning  in  Schoenberg  ([9]  Lemma  2 and  Theorem  3,  p.  260). 

Remark ■ The  result  of  Theorem  2 can  be  succintly  summarized  by  saying  that  the  map 

x -*•  Mx)  where  X(x)  = X^(x),  x e ( j— 1 , j I # j = l»...,n  is  a cn  1 map  of  (0,n)  onto 

(-<»,0)  . The  smoothness  of  this  map  follows  because  A (x;X)  is  a cardinal  L-spline. 

n 

Corollary  1.  If  T = T and  if  | Im  t^|  < m,  v = 0,  l,...,n,  then  for  any  complex  number  a 

the  roots  of  A (x;X|a  + T)  are  eaX. (x) , . . . ,eaX  (x)  where  {X  (x) }?  are  the  zeros  of 

n 1 n v 1 

An(xsX|T) . 

This  corollary  is  an  immediate  consequence  of  Theorem  3 and  follows  easily  on  using  the 


identity 


(a  + tQ,  a + t ,...,a  + tn]f(t)  = ItQf . . . ,tnl  f (t+a) 
and  the  representation  (4.4)  which  yields 

An(x;X|a  + T)  = ea(X  1)An(x;  e 3X|t)  . 

Theorem  4.  For  any  set  A of  n + 1 real  numbers  X.  < X.  <...<  X with  X - X.  < 2m  and 
‘ — 0 — 1 — — n n 0 

with  Xj  + X^  independent  of  j(j  = 0,  l,...,n),  the  polynomial  Hn(A;x|iA)  has  distinct 
roots  on  the  ray 


i(X  +X  )/2 

(6.6)  {-pe  0 |p  > 0}  . 

Moreover  II  (X  ,-x  I i A ) vanishes  for  some  X e U if  and  only  if  x = 0 and  n is  even  or 
n ...  ~ 

x = | and  n is  odd  and  in  each  of  these  cases  II(A;x  iA)  vanishes  only  for 

i(m+(X -+XJ/2) 

(6.7)  X = e 0 " 

Proof . The  first  part  of  the  theorem  follows  from  corollary  1 and  Theorem  3,  since  iA  = a + T, 

xQ+xn 

where  a = i — ^ — and 


(x  -x  ) <x  -xn  ,) 

r . On  1 n-1 

T = {l  j ' 1 2 


(X  -X„) 
n 0 


The  second  part  follows  from  the  observation  that  T = T = -T  which  implies,  because  of  (4.3) 
and  (4.6)  , that 

A (-i-; —1  It)  = 0 for  n odd,  A ( 0 ; —1 1 T)  = 0 for  n even 
n 2 n 

Corollary  2.  Given  any  m real  numbers  with  max|y  | < it  and  data  {y^}^  of  power 

v v 

growth  Y,  there  is  a unique  cardinal  L-spline  S(x)  of  power  growth  y corresponding  to  the 

1,1  2 2 

differential  operator  D II  (D  +u . ) , satisfying 

j=l  3 

S(v  + 6)  = yv,  v = 0,  +1,  +2,.... 
provided  B e (0,1) . When  B = 0,  this  result  fails. 

Proof . According  to  [6] , the  cardinal  L-spline  interpolation  problem  above  has  a unique  bound- 
► 

ed  solution  if  and  only  if 

n2m(^;S|iA)  + 0,  A = {iVjlj  = l,...,m}  U {0}  , 

for  A e U.  Since  A satisfies  the  conditions  of  Theorem  4,  the  corollary  is  seen  to  follow 
imnediately  from  it. 

Example  1.  In  Theorem  4,  set  A v = U+v)n,  i real,  |n|  < Then  An(x;A)  has  distinct 

zeros  on  the  ray 

+i  U + ^)n 

{-?e  ] C > 0} 

This  case  was  proved  by  Tzimbalario  [13) . 

Example  2 . If  we  take  A^  = vg,  with  nn  < 2tt  in  Theorem  4,  we  obtain  from  the  representa- 
tion (4.4)  with  x = 0 that 

(6.8)  An(0;A)  = [0,  in,  — ,inn)  — — , 

e -A 

has  n - 1 distinct  zeros  on  the  ray  {-c  expf^y^lc  > 0 . The  second  part  of  Theorem  4 im- 
plies that  when  n is  odd,  (6.8)  will  not  vanish  for  A = -1,  while  for  n even,  (6.8)  will 
vanish  for  A = -1  only  when  n * 0.  If  we  set  ein  = q,  then  for  n 4 0 an  easy  computa- 
tion shows  that  An(0;-1)  becomes  R^lq)  (except  for  a non-zero  factor)  where 


v=0  ' ' q +1 


Thus  we  get  the  result  that  Rn(q)  is  not  zero  on  the  arc  of  the  unit  circle  |z|  = 1 given 
by  | arg  q|  < ~ when  n is  odd  unless  a + 1.  When  n is  even,  Rn<q)  will  vanish  only 
when  n = 0,  i.e.,  q = 1.  In  [6]  it  was  shown  that  when  q is  real,  Rn<<?)  has  exactly 

n - 1 simple  positive  zeros  occuring  in  reciprocal  pairs.  For  the  imoortance  of  Rn(q)  in 
spline  interpolation  see  [2]  and  [6] . 

7.  Proof  of  Theorem  1,  According  to  Lemma  3 and  Theorem  4 (with  A replaced  by  aA) , the 
X.P  is  solvable  when  n is  even  and  t ^ 0 or  when  n is  odd  and  T 4 In  the  exceptional 

cases,  i.e.,  n even  and  T = 0 or  n odd  and  x — ~,  the  I.P  is  solvable  if  and  only  if 


(7.1) 


X„+X 
. 0 n 
i — 5 — a 

e e1”  + wV,  (v  = 0,  1, . . . ,k-l)  . 


X +X 
0 n 


This  implies  that  the  only  zero  of  II  (A;0  iaA)  on  the  unit  circle,  viz.  exp(i — - — a + iir)  is 

n * 

not  equal  to  a root  of  unity*  Since  oi  = u>  - e , it  follows  from  (7.1)  that 


(7.2) 


X ♦ x + k - 2v  * 0 (mod  2k),  v = 0,  l,...,k-l  . 


if  Aq  + X^  + k is  odd,  (7.2)  is  valid  so  that  the  I.P.  is  solvable.  This  completes  the  proof 
of  Theorem  1 . 

Proof  of  Theorem  2.  If  XQ  + Xn  + k is  even,  we  shall  show  that  (2.2)  fails  only  for 

v = (k+XQ+Xn)/2  = h (say) . Let  us  again  verify  what  we  saw  earlier  (see  the  proof  of  Lemma  2) 

that  if  v = X , (7.2)  is  equivalent  to  X -X  + k i 0 (mod  2k)  since  X + X 

1 n-r  r r n-r 

Aq  + An>  This  is  obviously  true  since  | Xn  f - Xf|  < k-1. 

If  Xr  < v < Af+1  for  some  r,  then  we  have 

X.  - X , +k<X„+X  +k-2v<X  -X  +k 

n-r-1  r+1  0 n n-r  r 

Since  X ’s  are  < k,  we  get 

v 

1<X„  + X +k-2v<  2k- 1 
0 n 

so  that  (7.2)  is  true  for  X < v < X 

r r+1 
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i 
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If  v > X , then 
n 

-k  < X + X -k  + 2<X  +X  + k - 2v  < X - X -2+k<k 
On  On  On 

Thus  X + X + k - 2v  = 0 (mod  2k)  if  and  only  if  X„  + X + k - 2v  = 0.  Thus  if  X + X + k 
On  ■'On  On 

is  even,  there  is  exactly  one  v for  which  inequality  (7.2)  fails.  From  Lemma  3 and  the  pre- 
ceding discussion  it  follows  that  (iji)  =0  if  and  only  if  j = h. 

If  S(z)  e S satisfies  (2.2),  then  by  Lemma  1, 

k-1 

S (z)  = l c.S.(z) 


so  that 


k-1  k-1 

l c S (i(iuV)  = l c.uV3S.  W = w , j = 0,  1, . . . ,k-l  . 


Therefore  c,S.(i|>)  = £ id  from  which  Theorem  2 easily  follows. 


Remark.  We  do  not  know  what  happens  if  + j '‘s  not  independent  of  j.  In  that  case 
the  above  method  does  not  seem  to  work. 

8.  Monosplines  and  quadrature  formulae. 

Consider  the  quadrature  formula 

k-1 

(8.1)  J f(z)dz  = £ c . f (ti)^ ) + Rf 

U j=0  3 

where  the  unit  circle  U is  described  counter-clockwise.  The  requirement  that  the  remainder 

X 

Rf  vanishes  when  f = z , v = 0,  l,...,n,  is  equivalent  to  the  relations 

X k-1  jX 

(8.2)  f z vdz  = 0 = £ c.w  V,  v = 0,  l,...,n 


We  assume  that  n + 1 < k in  which  case  there  are  k - n - 1 free  parameters  in  (8.1) . 

In  order  to  bring  out  the  connection  of  the  quadrature  formula  (8.1)  satisfying  (8.2) 
with  the  class  of  splines  S,  we  introduce  some  differential  operators.  Set 


i 


SB® 


V « D(* 


■ (A  -A  -1) 
V v-1 


f ) , v = 1,  2, . . . ,n 


£Qf  = D(z  uf) 


■(A  -A  -1) 


X f = z 

v 


Df,  XQf  = z Df,  v = 1,  2,...,n 


We  shall  prove 


A +1  n 

i £n  z " = n <x  -A.+l)  . 

n_1  0 j=0  " 3 


Lemma  5.  If 


- S(z)  , S(z)  e 5 , 


n ( a -a  +i) 

o n v 


and  if  f(z)  e C (U) , then 


/ f(z)dz  = l C.f(uj)  + (-l)n+1  / K ( z)  £*£*••  *X*f  dz  . 
U j=0  3 U 0 1 " 


Proof . It  is  easy  to  verify  that 

. * * k— 1 w3*1  -(A  -A  -1) 

/ K(z)  (XX  •••£  f)  (z)dz  = (-1)"  l f (X  •••Xk)(z)z  n n_i  f ' (z)dz  . 

0 0 1 n j=0  j n_1  0 

U) 

Since  H(z)  = z n n (Xn  • * *X0K)  (z)  has  discontinuities  at  z = u-1  (j  = 0,  l,...,n-l), 

the  coefficients  in  (8.4)  are  given  by  mu'1)  - H(w^). 

Formula  (8.4)  establishes  a correspondence  between  the  function  K(z)  which  we  will  call 

a A-monospline  and  quadrature  formulae  (8.1)  satisfying  (8.2) . 

A natural  problem  which  now  arises  is  to  find  a A-monospline  which  has  least  L -norm 

P 

1 < p < «.  We  seek  to  minimize 


|k||  p * / |K(z)  |P|dz| 


over  all  A-monosplines  of  order  n.  The  cuadrature  formula  which  corresponds  to  a solution 
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r* 


of  (8.5)  will  be  called  an  "optimal"  quadrature  formula. 

We  shall  prove 

Theorem  5.  Let  1 <_  n + 1 < k.  Among  all  quadrature  formulae  of  the  form  (8.1)  which  are 
exact  for  f e an  "optimal"  quadrature  formula  is 


(8.6) 


with 


, . k-1  (X  -n+1)  j 

/ f(z)dz  = ^ X,  l u n f(u3)  + Rf 

U k j=0 


r>4- 1 - * * 

Rf  = (-1)"  / K„(z)  (£0---JCnf(z))dz 


where  the  kernel  K4(z)  is  given  by 


K*(z)  = {zr  - AS  (z)},  r 5 (X  +1)  mod  k 

* n * r n 


H(r-X  ) 
0 


and  Xt  minimizes  the  integral 


/ | 1 - Xz  S (z) | d6  , 

o r 


ie 


over  all  X e I. 


Proof.  Let  K(z)  be  any  A -monospline . Set 


1 kr1  -£r  l 
K(z)  = i l id  “ K(a>  z) 


11=0 


Then  K(z)  is  a A-monospline  satisfying  the  functional  eouation 


Using  this  equation  the  theorem  easiiv  follows. 

Remark . When  p = 2 it  is  a straightforward  matter  to  compute  X^.  It  should  also  be  possible 
to  find  X#  when  p = 1,  ».  In  the  case  X^  = v considered  by  Schoenberg  (81,  Theorem  5 
leads  to  a solution  of  a question  raised  at  the  end  of  his  paper. 

9.  Conclusion. 

It  is  of  some  interest  to  point  out  that  Theorem  1 is  only  a sufficient  condition  for  the 

I.P  to  be  uniauely  solvable  in  the  class  S.  It  can  be  easily  proved  that  if  k is  even, 

i|»  “ 1,  Xg  « 0,  Xj  •=  1,  X2  = 4,  then  the  number  of  Lemma  3 given  by 

oo 

m _ y i 

r Um-m  (ks+r)  (ks+r-1)  (ks+r-4'  ' 

is  not  equal  to  zero.  This  is  easily  seen  on  rewriting 

OD  1 

Nr  * (ks+r) (ks+r-1) (ks+r-4) 
s-0 

00 

_ y i j . i 

(ks+k-r+4) (ks+k-r+l) (ks+k-r)  1 2 

s«0 

A comparison  of  the  terms  in  1^  and  I2  shows  that  for  2r  <_  k+2 , 1^  > I2  each  term  in  1^ 
being  greater  than  the  corresponding  term  in  Ij.  Similarly  for  2r  >_  k+2,  1^  < I2>  Thus  the 
only  case  to  be  checked  is  2r  = k+3  which  cannot  occur  since  k is  even. 

Thus  although  XQ  + X2  ^ 2X^,  the  criterion  (4.13)  of  Lemma  3 applies  and  the  I.P  is 
uniquely  solvable.  It  would  be  interesting  to  obtain  a proof  of  Theorem  1 independent  of  the 
theory  of  Cardinal  L-splines  and  to  see  if  the  conditions  on  A can  !e  further  relaxed. 
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